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-Abstract- 

The Continuous Skolem Problem asks whether a real-valued function satisfying a linear differen¬ 
tial equation has a zero in a given interval of real numbers. This is a fundamental reachability 
problem for continuous linear dynamical systems, such as linear hybrid automata and continuous¬ 
time Markov chains. Decidability of the problem is currently open—indeed decidability is open 
even for the sub-problem in which a zero is sought in a bounded interval. In this paper we show 
decidability of the bounded problem subject to Schanuel’s Conjecture, a unifying conjecture in 
transcendental number theory. We furthermore analyse the unbounded problem in terms of the 
frequencies of the differential equation, that is, the imaginary parts of the characteristic roots. 
We show that the unbounded problem can be reduced to the bounded problem if there is at most 
one rationally linearly independent frequency, or if there are two rationally linearly independent 
frequencies and all characteristic roots are simple. We complete the picture by showing that de¬ 
cidability of the unbounded problem in the case of two (or more) rationally linearly independent 
frequencies would entail a major new effectiveness result in Diophantine approximation, namely 
computability of the Diophantine-approximation types of all real algebraic numbers. 
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[T] Introduction 

The Continuous Skolem Problem is a fundamental decision problem concerning reachability in 
continuous-time linear dynamical systems. The problem asks whether a real-valued function 
satisfying an ordinary linear differential equation has a zero in a given interval of real numbers. 
More precisely, an instance of the problem comprises an interval I C K>q with rational 
endpoints, an ordinary differential equation 

/(") + + ... + ao/= 0 (1) 

whose coefficients are real algebraic, together with initial conditions /(O),..., 
that are also real algebraic numbers. Writing / : M>o —>■ K for the unique solution of the 
differential equation subject to the initial conditions, the question is whether there exists 
t £ I such that f{t) = 0. Decidability of this problem is currently open. Decidability of 
the sub-problem in which the interval I is bounded, called the Bounded Continuous Skolem 
Problem, is also open [H Open Problem 17]. 

The nomenclature Continuous Skolem Problem is based on an analogy with the Skolem 
Problem for linear recurrence sequences, which asks whether a given linear recurrence 
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sequence has a zero term jl2j . Whether the latter problem is decidable is an outstanding 
question in number theory and theoretical computer science; see, e.g., the exposition of 
Tao Uni Section 3.9]. 

The continuous dynamics of linear hybrid automata and the evolution of continuous¬ 
time Markov chains, amongst many other examples, are determined by linear differential 
equations of the form x'{t) = Ax(t), where x{t) G K" and A is an n x n matrix of real 
numbers PQ. A basic reachability question in this context is whether, starting from an initial 
state a;(0), the system reaches a given hyperplane {y G K" : u^y = 0} with normal vector 
u G M". For example, one can ask whether the continuous flow of a hybrid automaton 
leads to a particular transition guard being satisfied or an invariant being violated. Now 
the function f{t) — v?"x{t) satisfies a linear differential equation of the form Q, and it 
turns out that the hyperplane reachability problem is inter-reducible with the Continuous 
Skolem Problem (see [H Theorem 6] for further details). Moreover, under this reduction 
the Bounded Continuous Skolem Problem corresponds to a time-bounded version of the 
hyperplane reachability problem. 

The characteristic polynomial of the differential equation 0 is 

X{x) := x” -I- ao . 

Let Ai,...,Am be the distinct roots of y. Any solution of 0 has the form f{t) = 
where the Pj are polynomials with algebraic coefficients that are determined 
by the initial conditions of the differential equation. We call a function / in this form 
an exponential polynomial. If the roots of y are all simple then / can be written as an 
exponential polynomial in which the polynomials Pj are all constant. 

The Continuous Skolem Problem can equivalently be formulated in terms of whether an 
exponential polynomial has a zero in a given interval of reals. If the characteristic roots 
have the form Xj = rj + iujj, where rj,ujj G K, then we can also write f(t) in the form 
f{t) = + Q 2 ,j{t) cos{ujjt)), where the polynomials Qij,Q 2 ,j have 

real algebraic coefficients. We call wi,... ,u!m the frequencies of /. 

Our first result is to show decidability of the Bounded Continuous Skolem Problem subject 
to Schanuel’s Conjecture, a unifying conjecture in transcendental number theory that plays a 
key role in the study of the exponential function on both real and complex numbers 
Intuitively, decidability of the Bounded Continuous Skolem Problem is non-trivial because 
an exponential polynomial can approach 0 tangentially. Assuming Schanuel’s Conjecture, 
we show that any exponential polynomial admits a factorisation such that the zeros of each 
factor can be detected using finite-precision numerical computations. Our method, however, 
does not enable us to bound the precision required to find zeros, so we do not obtain a 
complexity bound for the procedure. 

A celebrated paper of MacIntyre and Wilkie m obtains decidability of the first-order the¬ 
ory of Kexp = 0,1, <, •, -f, exp) assuming Schanuel’s Conjecture over M. The proof of [171 

Theorem 3.1] mentions an unpublished result of MacIntyre and Wilkie that generalises [IH] 
to obtain decidability when Kexp is augmented with the restricted functions sin [[o, 27 r] and 
cos [[o, 27 r]i tills time assuming Schanuel’s Conjecture over C. This result immediately implies 
(conditional) decidability of the Bounded Continuous Skolem Problem. However, decidability 
of latter problem is simpler and, as we show below, can be established more directly. 

In the unbounded case we analyse exponential polynomials in terms of the number of 
rationally linearly independent frequencies. We show that the unbounded problem can be 
reduced to the bounded problem if there is at most one rationally linearly independent 
frequency, or if there are two rationally linearly independent frequencies and all characteristic 
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roots are simple. These two reductions are unconditional and rely on the cell decomposition 
theorem for semi-algebraic sets [3] and Baker’s Theorem on linear forms in logarithms of 
algebraic numbers [2]. 

We complete the picture by showing that decidability of the unbounded problem in the 
case of two (or more) rationally linearly independent frequencies would entail a major new 
effectiveness result in Diophantine approximation—namely computability of the Diophantine- 
approximation types of all real algebraic numbers. As we discuss in Appendix]^ currently 
essentially nothing is known about Diophantine-approximation types of algebraic numbers of 
degree three or higher, and they are the subject of several longstanding open problems. 

The question of deciding whether an exponential polynomial / has infinitely many zeros 
is investigated in There the problem is shown to be decidable if / satisfies a differential 
equation of order at most 7. This result does not rely on Schanuel’s Conjecture. It is also 
shown in [7] that, analogously with the Continuous Skolem Problem, decidability of the 
Infinite Zeros Problem in the general case would entail significant new effectiveness results in 
Diophantine approximation. 

Mathematical Background 

2.1 Zero Finding 

Let / : [a, 6] —>■ K be a function defined on a closed interval of reals with endpoints a,b G Q. 
Suppose the following two conditions hold: (i) there exists M > 0 such that / is M-Lipschitz, 
i.e., |/(s) — f{t)\ < M\s — t\ for all s,t G [a, &]; (ii) given t G [a, 6] C Q and positive error 
bound £ S Q, we can compute q G Q such that \f{t) — q\ < e. Then given a positive 
rational number S we can compute piecewise linear functions : [a, 6] —>■ M such that 

fi{t) < f{t) < f^i^) s-iid f/{t) — f^{t) < S for all t G [a, 6]. We do this as follows: 

1. Pick N G N such that ^ and consider sample points Sj := a + 

j = 0,..., A^, dividing the interval [a, b] into N sub-intervals, each of length at most 

2 . For each sample point Sj compute qj G Q such that \qj — f{sj)\ < define 
f^{sj) = qj + and extend and linearly between sample points. 

Note that the Lipschitz condition on / ensures that < f < ■ 

Now suppose that / satisfies the following additional conditions: (iii) /(a) ^ 0, f{b) ^ 0; 
(iv) for any t G (a, b) such that f{t) = 0, f'{t) exists and is non-zero, i.e., / has no tangential 
zeros. Then we can decide the existence of a zero of / by computing upper and lower 
approximations and for successively smaller values of 5. If (t) < 0 for all t or 
f^{t) > 0 for all t then we conclude that / has no zero on [a, b]; if (s) < 0 and fs{t) > 0 
for some s, t then we conclude that / has a zero; otherwise we proceed to a smaller value 
of S. This procedure terminates since by (iii) and (iv) either / has a zero in [a, b] or it is 
bounded away from zero. 

2.2 Number-Theoretic Algorithms 

For the purposes of establishing decidability, we can assume that an instance of the Continuous 
Skolem Problem is a real-valued exponential polynomial f{t) = where 

Ai,..., Am and the coefficients of the polynomials Pi,..., Pm are algebraic, see [D Theorem 

6 ]. 

For computational purposes we represent an algebraic number a by a polynomial P with 
rational coefficients such that P{a) = 0, together with a numerical approximation p + qi, 
where p,q G Q, of sufficient accuracy to distinguish a from the other roots of P [HI Section 
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4.2.1]. Given this representation we can obtain approximations of a to within an arbitrarily 
small additive error. 

Let K be the extension field of Q generated by Ai,..., Xm and the coefficients of the 
polynomials Pi,..., Pm- Note that K is closed under complex conjugation. We can compute 
a primitive element of K, that is, an algebraic number 6 such that K = Q(ff), together with 
a representation of each characteristic root Xj as a polynomial in 6 with rational coefficients 
(see |H1 Section 4.5]). From the representation of Ai,...,Am as elements of Q{0), it is 
straightforward to determine maximal Q-linearly independent subsets of {Re(Aj) : 1 < j < m} 
and {Im(Aj) : 1 < j < m} (see [HI Section 1]). 

Let log denote the branch of the complex logarithm defined by log(re*®) = log(r) + iO 
for a positive real number r and 0 < 0 < 2tt. Recall that one can compute log z and to 
within arbitrarily small additive error given a sufficiently precise approximation of z |B]. 

2.3 Laurent Polynomials 

Let RT be a sub-field of C that has finite dimension over Q and is closed under complex 
conjugation. Fix non-negative integers r and s, and consider a single variable x and tuples 
of variables y = (yi,..., y^) and z = (zi,..., Zg}- Consider the ring of Laurent polynomials 

n ■= R:[a;,yi,y]"\...,yr,y;:\zi,zf\...,zs,z7^] , 

which can be seen as a localisatioi|^ of the polynomial ring A := K[x^ yi,..., y^, zi,..., Zs] 
in the multiplicative set generated by the set of variables {yi,..., y^} U {zi,..., Zs}. The 
multiplicative units of TZ are the non-zero monomials in variables yi,..., y^ and zi,..., Zg. As 
the localisation of a unique factorisation domain, TZ is itself a unique factorisation domain [21 
Theorem 10.3.7]. From the proof of this fact it moreover easily follows that TZ inherits from 
A computability of factorisation into irreducibles (e.g., using the algorithm of m)- 

We extend the operation of complex conjugation to a ring automorphism of TZ as follows. 
Given a polynomial 


n 



where oi,..., a„ € AT, define its conjugate to be 


n 



This definition is motivated by thinking of the variables x and yi,..., y^ as real-valued and 
the variables Zi,..., Zg as taking values in the unit circle in the complex plane. 

We will need the following proposition characterising those polynomials in P S 7?. such 
that P and P are associates, i.e., such that P is equal to the product of P by a monomial. 
Here we use pointwise notation for exponentiation: given a tuple of integers u = (ui,..., Mg), 
we write z“ for the monomial z“^ ... Zg “. 

► Proposition 1. Let P G TZ he such that P = z^P for it G Then either (i) P has 
the form P = z'^Q for some Q G TZ with Q = Q, or (ii) there exists Q G TZ such that 
P — Q + z'^Q and P does not divide Q in TZ. 


^ Recall that the localisation of a commutative ring W in a multiplicatively closed subset S such that 
Ou ^ S is the ring of formal fractions Us = {c/s ■ a GU,s G S}, with addition and multiplication 
defined as usual. 
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Proof. Consider a monomial M such that z'^M = M. Then M has a real coefficient and the 
exponent ii; of 2 : in M satisfies 2w = u. Thus if z'^M = M for every monomial M appearing 
in P then P has the form Qz^, where 2w = u and Q is a polynomial in the variables x and 
y with real coefficients. In particular Q = Q, and statement (i) of the proposition applies. 

Suppose now that z'^ M ^ M for some monomial M appearing in P. Then the map 
sending M to z'^M induces a permutation of order 2 on the monomials on P. Thus we may 
write P = where n = k + 2^ for some k > 0 and £ > 1 such that z'^Mj = Mj for 

^ < j < k and z'^Mj = Mj+£ for k+1 < j < i. Then, writing Q := 5 J2'j=i 
we have P = Q + z'^Q. 

The set of monomials appearing in Q is a proper subset of the set of monomials appearing 
in P (up to constant coefficients) and so Q cannot be a constant multiple of P. It also follows 
that for each variable a G {a;, yj, Zk ■ ^ < j < r,l < k < s}, the maximum degree of tr in P 
is at least its maximum degree in Q, and likewise for a~^. Thus Q cannot be a multiple of 
P by a non-constant polynomial either. We conclude that P does not divide Q. ◄ 


2.4 Transcendence Theory 

We will use transcendence theory in our analysis of both the bounded and unbounded variants 
of the Continuous Skolem Problem. In the unbounded case we will use the following classical 
result. 


► Theorem 2 (Gelfond-Schneider). Let a,b be algebraic numbers not equal to 0 or 1. Then 
for any branch of the logarithm function, is either rational or transcendental. 

In fact we will make use of the following corollary, which is obtained by applying Theorem 
to the algebraic numbers a = and b = 


► Corollary 3. Let ai yf fSi, a 2 yf (^2 all lie in [ 0 , 7 r] and suppose t/iat cos(ai), cos(a 2 ), cos(/3i) 

and cos(/32 ) are algebraic. Then is either rational or transcendental. 

Our results in the bounded case depend on Schanuel’s conjecture, a unifying conjecture 
in transcendental number theory which, if true, greatly generalises many of the cent¬ 
ral results in the field (including the Gelfond-Schneider Theorem, above). Recall that a 
transcendence basis of a field extension L/K is a subset S C L such that S is algebraically 
independent over K and L is algebraic over K{S). All transcendence bases of L/K have the 
same cardinality, which is called the transcendence degree of the extension. 

► Conjecture 4 (Schanuel’s Conjecture |15j). Let ai,...,a„ be complex numbers that are 
linearly independent over Q. Then the field Q(ai,..., a„, e“L ..., e“") has transcendence 
degree at least n over Q. 

A special case of Schanuel’s conjecture, that is known to hold unconditionally, is the 
Lindemann-Weierstrass Theorem m- if oi,..., a„ are algebraic numbers that are linearly 
independent over Q, then e“L ■ • ■ i 6 “" are algebraically independent. 

We apply Schanuel’s conjecture via the following proposition. 

► Proposition 5. Let {ai,...,ar} and {&i,...,5s} be Q-linearly independent sets of real 
algebraic numbers. Furthermore, let P,Q G TZ he two polynomials that have algebraic 
coefficients and are coprime in TZ. Then the equations 


P(t,e“i‘,.. 

a^t Ab-i t 

.,e " ,e W-- 


Q(t,e“i‘,.. 

art ibit 



have no non-zero common solution t G IR. 


( 2 ) 

( 3 ) 
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Proof. Consider a solution t ^ 0 of Equations Q and (|^. By passing to suitable associates, 
we may assume without loss of generality that P and Q lie in A, i.e., that all variables in P 
and Q appear with non-negative exponent. Moreover, since P and Q are coprime in TZ, their 
greatest common divisor i? in M is a monomial. In particular. 


Thus, dividing P and Q by R, we may assume that P and Q are coprime in A and that 
Equations Q and ([^ still hold. 

By Schanuel’s conjecture, the extension 


Q(ait,..., art, ibit,ibst, 6“^*,..., e 


aj^t AbTt 

e S • ■ ■ ,e 


has transcendence degree at least r -|- s. Since oi,...,and hi,... ,bs are algebraic over Q, 
writing 

it follows that the extension Q(S')/Q also has transcendence degree at least r -|- s. 

From Equations and (§ we can regard S as specifying a common root of P and Q. 
Pick some variable a S {x, Vj, Zj '. 1 < i < r,l < j < s\ that has positive degree in P. Then 
the component of S corresponding to cr is algebraic over the remaining components of S. We 
claim that the remaining components of S are algebraically dependent and thus S comprises 
at most r -I- s — 1 algebraically independent elements, contradicting Schanuel’s conjecture. 
The claim clearly holds if a does not appear in Q. On the other hand, if a has positive 
degree in Q then, since P and Q are coprime in A, the multivariate resultant ReSo-(T’, <5) is 
a non-zero polynomial in the set of variables {x, yj, Zj : 1 < i < r,l < j < s} \ {cr} which 
has a root at S (see, e.g., IM Page 163]). Thus the claim also holds in this case. In either 
case we obtain a contradiction to Schanuel’s conjecture and we conclude that Equations ([^ 
and ^ have no non-zero solution t G K. ◄ 


Decidability of the Bounded Continuous Skolem Problem 

Suppose that {ai,...,ar} and {ibi,... ,ibs} are Q-linearly independent sets of real and 
imaginary numbers respectively. Let the ring of Laurent polynomials TZ be as in Section |2.3| 
and consider the exponential polynomial 

/(t) = P(t,e“^‘,...,e“’-‘,e*'^*,...,e*''=‘), (4) 

where P G TZ is irreducible. We say that / is a Type-1 exponential polynomial if P and P 
are not associates in TZ, we say that / is Type-2 ii P = aP for some a G C, and we say that 
/ is Type-3 if P = t/P for some non-constant unit U G TZ. 

► Example 6. The simplest example of a Type-3 exponential polynomial is g{t) = 1 -|- 
Here g{t) = P(e®‘), where P{z) = 1 -I- x is an irreducible polynomial that is associated with 
its conjugate P{z) = 1 -|- Note that the exponential polynomial f{t) — 2 + 2cos(t), 
which has infinitely many tangential zeros, factors as the product of two type-3 exponential 
polynomials f{t) = g{t)g{t). 

In the case of a Type-2 exponential polynomial P = aP it is clear that we must have 
jaj = 1. Moreover, by replacing P by /3P, where /3^ = a, we may assume without loss of 
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generality that P = P. Similarly, in the case of a Type-3 exponential polynomial, we can 
assume without loss of generality that P = z'^P for some non-zero vector tt G Z'*. 

Now consider an arbitrary exponential polynomial /(t) := Assume that 

the coefficient field K of TZ contains the coefficients of Pi,..., Pm- Let {oi,..., 0 ,.} be a 
basis of the Q-vector space spanned by {Re(Aj) '■ 1 < j < m} and let {&i,..., bs} be a basis 
of the the Q-vector space spanned by {Im(Aj) '■ 1 < j < m}. Without loss of generality we 
may assume that each characteristic root A is an integer linear combination of oi,..., 
and ibi,. ■ ■ jibs- Then is a product of positive and negative powers of ..., e“''‘ and 
..., and hence there is a Laurent polynomial P G TZ such that 


/(t) = P(t,e“i‘ 


^a-rt ^ibit 




)• 


( 5 ) 


Since P can be written as a product of irreducible factors, it follows that / can be written 
as product of Type-1, Type-2, and Type-3 exponential polynomials, and moreover this 
factorisation can be computed from /. Thus it suffices to show how to decide the existence 
of zeros of these three special forms of exponential polynomial. We will handle all three cases 
using Schanuel’s conjecture. 

Writing the exponential polynomial /(<) in (j^ in the form f{t) = : if 

follows from the irreducibility of P that the polynomials Qi,, Qm have no common root. 
But then by the Lindemann-Weierstrass Theorem any zero of / must be transcendental 
(see m Theorem 8]). 

► Theorem 7. The Bounded Continuous Skolem Problem is decidable subject to Schanuel’s 
conjecture. 

Proof. Consider an exponential polynomial 

/(t) = P(t,e“^‘,...,e“’-‘,P'’^‘,...,P''=*), (6) 

where P G TZ is irreducible. Suppose that {ai,... ,ar} and {ibi,... ,ibs} are Q-linearly 
independent sets of, respectively, real and imaginary numbers lying in the coefficient field 
K of TZ. We show how to decide whether / has a zero in a bounded interval I C K>o, 
considering separately the case of Type-1, Type-2, and Type-3 exponential polynomials. 


Case (i): / is a type-1 exponential polynomial 

Note that P and P are coprime in TZ since, by assumption, they are both irreducible and 
are not associates. We claim that in this case the equation fit) = 0 has no solution t gM.. 
Indeed f{t) = 0 implies 

and the non-existence of a zero of / follows immediately from Proposition!^ 


Case (ii): / is a type-2 exponential polynomial 

In this case we have P = P and so / is real-valued. Our aim is to use the procedure of 
Section 2.1 to determine whether or not / has a zero in [c, d], where c,d G Q. To this 
end, notice first that f{c),f{d) ^ 0 since any root of / must be transcendental. Moreover, 
since f is bounded on [c, d], / is Lipschitz on [c, d]. It remains to verify that the equations 
fit) = 0, fit) = 0 have no common solution t G [c, d]. 
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We can write f'{t) in the form 


where Q is the polynomial 


Q 


dx 


i=i i=i 


dP 

dzj ■ 


We claim that P and Q are coprime in P. Indeed, since P is irreducible, P and Q can only 
fail to be coprime if P divides Q. 

If P has strictly positive degree fc in a; then Q has degree fc — 1 in a; and thus P cannot 
divide Q. (Recall that all polynomials in TZ have non-negative degree in the variable x.) 
On the other hand, if P has degree 0 in a; then Q is obtained from P by multiplying each 
monomial y'^z'" appearing in P by the complex-number constant 

Moreover, by the assumption of linear independence of {oi,..., o^} and {6i,..., &s}, each 
monomial in P is multiplied by a different constant. Since P is not a unit, it has at least two 
different monomials and so P is not a constant multiple of Q. Furthermore, for each variable 
a G {yj,yj^ : 1 < j < r} U {zj,z~^ : 1 < j < s}, its degree in P is equal to its degree in Q. 
Thus P cannot be a multiple of Q by a non-constant polynomial either. 

We conclude that P does not divide Q and hence P and Q are coprime. It now follows 
from Propositionthat the equations f{t) = f'{t) = 0 have no solution < G K. 


Case (iii): / is a type-3 exponential polynomial 

Suppose that / is a Type-3 exponential polynomial. Then in ^ we have that P — z'^P 
for some non-zero vector m G Z'*. By Proposition we can write P = Q + z'^Q for some 
polynomial Q € TZ that is coprime with P. 

Now define 


gi{t) ■■= Q{t,e 


ait 




and g2{t) := gi{t), so that f{t) = gi{t) + g2{t) for all t. 

We show that 52 (^) 7^ 0 for all t G K. Indeed if g2{t) = 0 for some t then we also have 
gi{t) = 0 and hence f{t) = 0. For such a t it follows that 

But P and Q are coprime and so these two equations cannot both hold by Proposition!^ Not 
only do we have 32(^) 7^ 0 for all t G K, but, applying the sampling procedure in Section [ 2 T] 
to | 52 (f)P (which is a differentiable function) we can compute a strictly positive lower bound 
on 1 ^ 2 (^)I over the interval [c,d]. 

Since 32 (i) 7 ^ 0 for all t G M we may define the function h : [c, d] —> K by 


h{t) := -K + ilog 



Notice that h{t) = 0 if and only if f{t) = 0. Our aim is to use the procedure of Section 2.1 


to decide the existence of a zero of h in the interval [c, d], and thus decide whether / has a 
zero in [c, d]. 
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Let t e (c, d) be such that h{t) = 0. Then gi{t) = —g 2 {t) and so = —1 does not lie 
on the branch cut of the logarithm function. It follows that h is differentiable at t and 


h’it) = 0 


92{t) 9'i{t)g2{t) - g2{t)gi{t) ^ 

9i{t) 52 

iff 9 'i(t)g 2 {t) - g 2 {t)gi{t) = 0 (since | 5 i(t)| = | 52 (t)| 0 ) 

iff 9'i{t)92{t) + 52 W 52 (t) = 0 (since 51 (t) = - 52 (t)) 

iff 5i(^) + 52W = 0 

iff f'{t) = 0. 


Thus h{t) = h'{t) = 0 implies f{t) = f'{t) = 0. But the proof in Case (ii) shows that 
f{t) = f'{t) = 0 is impossible. (Nothing in that argument hinges on / being real-valued.) 
Thus h has no tangential zeros in (c, d). 

We cannot directly use the procedure in Section 2.1 to decide whether h has a zero in 
[c, d] since h is not necessarily continuous: its value can jump from —tt to tt (or vice versa) 
due to the branch cut of the logarithm along the positive real axis. However, due to the 
strictly positive lower bound on | 52 (t)|, the function \h\ is Lipschitz on [c,d]. Thus, applying 
the sampling procedure in Section |2.1| for computing lower and upper bounds of Lipschitz 
functions we can compute a set E C [c, d] such that if is a finite union of intervals with 
rational endpoints, |/(t)| < ^ for t G if, and \f{t)\ > | for t ^ if. In particular, E contains 
all zeros of / in [c, d] and / is Lipschitz on E. Thus we can apply the zero-finding procedure 
from Section |2.1| to the restriction h \ E and thereby decide whether or not h has a zero on 
[c,d]. 


[~^ The Unbounded Case 

In this section we consider the unbounded case of the Continuous Skolem Problem. For our 
analysis it is convenient to present exponential polynomials in the form 

n 

fit) = i^hjit) cos{u}jt) + P 2 ,j{t) sin{uj jt)) , (7) 

1=1 

where rj,ujj are real algebraic numbers and Pij,P 2 ,j are polynomials with real algebraic 
coefficients for j = 1,..., n. Our aim is to classify the difficulty of the problem in terms of 
the number of rationally linear independent frequencies wi,... ,a;„. 

Recall that in Section we have shown the bounded problem to be decidable subject to 
Schanuel’s Conjecture. In Appendix we give a reduction of the unbounded problem to the 
bounded problem in case the set of frequencies spans a one-dimensional vector space over Q. 
In the present section we give a reduction of the unbounded problem to the bounded problem 
in case the set of frequencies spans a two-dimensional vector space over Q and the polynomials 
Pij and P 2 J are all constant. (This last condition is equivalent to the assumption that f{t) 
is simple.) The argument in the two-dimensional case is a more sophisticated version of that 
in the one-dimensional case, although the result is not more general due the assumption of 
simplicity. 

In Appendix we present a family of instances showing that obtaining decidability of 
the unbounded problem in the two-dimensional case without the assumption of simplicity 
would require much finer Diophantine-approximation bounds than are currently known. 
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4.1 Background on Semi-Algebraic Sets 

A subset of K" is semi-algebraic if it is defined by a Boolean combination of constraints 
of the form P{xi ,..., Xn) > 0, where P is a polynomial with real algebraic coefficients. A 
partial function / : K” —)■ M is semi-algebraic if its graph is a semi-algebraic subset of 
The Tarski-Seidenberg Theorem O Section 1] states that the semi-algebraic sets are closed 
under projection and are therefore precisely the first-order definable sets over the structure 
(M, <, -|-, •, 0 , 1 ). 

Let {ii,... ,in) be a sequence of zeros and ones of length n > 1. An (ii,..., in)-ceU is a 
subset of K", defined by induction on n as follows: 

(i) A (O)-cell is a singleton subset of K and a (l)-cell is an open interval (a, b) C M. 

(ii) Let X C K” be a (A,... ,i„)-cell and / : A —>■ M a continuous semi-algebraic function. 

Then {{x,f{x)) S : x G X} is a (ii,..., 0)-cell, while {{x,y) G : x G 

X Ay < f{x)} and {{x, y) G : x G X Ay > f{x)} are both (ii,..., l)-cells. 

(iii) Let X C M” be a (A,..., i„)-cell and /, g : A —> K continuous semi-algebraic functions 

such that f{x) < g{x) for all x G X. Then {{x,y) G : f{x) < y < g{x)} is a 

(ii,.. .,in, l)-cell. 

A cell in M" is a (ii,...,*„)-cell for some (necessarily unique) sequence (ii,..., f„). 

A fundamental result about semi-algebraic sets, that we will use below, is the Cell- 
Decomposition Theorem [3]; given a semi-algebraic set E C M" one can compute a partition 
of P as a disjoint union of cells P = Ci U ... U Cm- 

We will also need the following result, proved in Appendix [B| 

► Lemma 8 . Let D C M" be a semi-algebraic set, g : P —>■ K a bounded semi-algebraic 
function, and ri,..., r„ real algebraic numbers. Define S = {t G M>o : (e'-i‘,...,e^"‘) GP}. 
Then 

(i) It is decidable whether or not S is bounded. If S is bounded then we can compute Pq G N 
such that S C [0, To] o-nd if S is unbounded then we can compute Tq G N such that 
(To,oo) C 5. 

(ii) If S is unbounded then the limit g* = limt_>oo ..., e’'"*) exists, is an algebraic 

number, and there are effective constants Ti, e > 0 such that \g(e '^^*,..., e’'"*) —g*\< 

for all t > Ti. 

4.2 Two Linearly Independent Frequencies 

The following lemma, which is a reformulation of |H Lemma 13], plays an instrumental role 
in this section. The lemma itself relies on a powerful quantitative result in transcendence 
theory—Baker’s Theorem on linear forms in logarithms of algebraic numbers |2]. 

► Lemma 9. Let 6 i ,&2 be real algebraic numbers, linearly independent over Q. Furthermore, 
let ipi,(p 2 be real numbers such that and are algebraic. Then there exist effectively 
computable constants N,T > 0 such that for all t > T and all ki,k 2 G Z, at least one of 
\bit — ipi — 2 fci 7 r| > 1/t^ and \b 2 t — <P 2 — ‘ 2 k 2 'x\ > 1/t^ holds. 

The main result of the section is the following. 

► Theorem 10. Let f(t) = (®i.i cos(wj<) -\- a 2 j sin{ujjt)) be an exponential poly¬ 

nomial where rj,aij, a 2 j,i 0 j are real algebraic numbers and the Q-span of {wi,..., w„} has 
dimension two as a Q-vector space. Then we can decide whether or not {t G R>o : f{t) = 0} is 
bounded and, if bounded, we can compute an integer T such that {t G K>o : fft) = 0} C [0, T]. 
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Proof. Let 6 i ,&2 be real algebraic numbers, linearly independent over Q, such that ojj is 
an integer linear combination of &i and 62 for j = 1, ... ,n. For each n G Z, sin(n&it) and 
cos{nbit) can be written as polynomials in sin( 6 it) and cos( 6 it) with integer coefficients, and 
similarly for 62 - It follows that we can write / in the form 

f{t) = , e’’"*,cos( 6 if), sin( 6 it),cos( 62 t),sin( 62 t)) 

for some polynomial Q with real algebraic coefficients that is computable from /. 

Write i?++ = {t > 0 : sin( 6 it) > 0 Asin(& 2 t) > 0}, R+- = {t > 0 : sin( 6 it) > 0 Asin( 62 t) < 

0}, and likewise define R _ R _for the two remaining sign conditions on sin( 6 it) and 

sin( 62 t). We show how to decide boundedness of {t G i?++ : f(t) = 0}. (The cases for 

i?_i_, R_and R_follow mutatis mutandis.) The idea is to compute a partition of 

{t G R++ : f{t) = 0} into components Zi,..., and to separately decide boundedness of 
each component Zj. 

Define a semi-algebraic set 

E = {{u,xi,X 2 ) G : 3yi, j /2 > 0 {xf + yf = xl + yj = I A Q{u,xi,yi,X 2 ,y 2 ) = O) } • 

Then for t G R++ we have f{t) = 0 if and only if (e’’‘,cos( 6 it),cos( 62 t)) G E, where r = 
(ri,..., Tn). Now consider a cell decomposition E = CiLi.. .UCm for cells Ci,..., Cm C 
and define 

Zj = {t G R++ : (e’’‘,cos(&if),cos( 62 t)) G Cj} , j = 1,... ,m, ( 8 ) 

Then {t G R++ '■ f{t) = 0} = Zi U ... U Zm- 

Now fix j G m}. We show how to decide boundedness of Zj. To this end, write 

Dj C K" for the projection of the corresponding cell Cj C on the first n coordinates. 

First suppose that {t G K : e’’* G Dj} is bounded. Then by Lemmawe can compute an 
upper bound T of this set. But Zj C {t G K>o : e'’* G Dj} and so Zj C [0,r]. 

On the other hand, suppose that {t G M : e’’* G Dj} is unbounded. Then, by Lemma 
this set contains an unbounded interval (T, 00 ) for some T G N. Write I = [—1,1] and define 
functions 51 , 52 , hi, h 2 : Dj —>■ K by 

gi{u) = \TE{x Gl ■.3y{u,x,y) gCj} 52 (m) = infjy G / : 3a: (m, a;, 5 ) G (9) 

/ii(m) = sup{a: G / : 35 (m, a:, y) G Cj} h 2 (M) = sup{y G / : 3a: (t 6 , a;, y) G Cj} (10) 

These functions are all semi-algebraic by quantifier elimination. Hence by Lemma the 
limits 5 * = limt_,,oo 5 i(e’^*) and h* = exist for i = 1,2 and are algebraic 

numbers. Clearly we have 5 } < h} and We now consider three cases according to 

the strictness of these inequalities. 

4.2.1 Case I: g} = h} and = h^. 

We show that Zj is bounded and that we can compute T 2 such that Zj C [ 0 ,r 2 ]. 

By Lemmathere exist Ti,e > 0 such that for all t > Ti and i = 1,2, 

l 5 .(e"‘) - 5*1 < e-^‘ and |h,(e"‘) - h*! < . ( 11 ) 

Then for t G R++ such that t > Ti we have 
t G Zj (e^\cos(bit),cos(b 2 t)) G Cj (by (|^) 

=> 5 i(e’’*) < cos( 6 it) < hi(e'’*) and 52 ( 6 ^*) < cos(b 2 t) < h 2 (e’'*) (by (|^(|To|)) 
=> |cos( 6 it) — 5 }| < and |cos( 62 t) — 52 I < (by pT] )) ( 12 ) 
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Write 5* = cos(v?i) and = cos((/32) for some </3i,<p2 G [0,7r]. Since |cos((/7i + x) — 
cos((/3i)| > x ^/3 for all x sufficiently small (by a Taylor expansion), the inequality ( 12 1 implies 
that for some fci, ^2 G Z, 


\bit — ipi — 2 ki'K\ < 3 e and \b2t — ip2 — 2fc27r| < 3 e . 


( 13 ) 


Combining the upper bounds in ( 13 ) with the polynomial lower bounds \bit — ipi — 2fci7r| > 
1 /t^ and |62t — 1^2 — 2A:27r| > 1 /t^ from Lemmawe obtain an effective bound T2 for which 
t G Zj implies t < T2. 


4.2.2 Case II: < h^. 

In this case we show that Zj is unbounded. The geometric intuition is as follows. We imagine 
a particle in the plane whose position at time t is (cos( 6 it),cos( 62 t)), together with a “moving 
target” whose extent at time t is Tt = {(x, y) : (e'’‘,x,y) G Cj}. Below we essentially argue 
that such a particle is bound to hit Tj at some time t since its orbit is dense in [— 1 ,+ 1 ]^ 
and Ft has positive dimension in the limit. 

Proceeding formally, first notice that Cj cannot be a (..., 0, l)-cell or a (..., 0,0)-cell, 
for then we would have gi{u) = hi{u) for all u G Dj and hence y* = h\. Thus Cj must 
either be a (..., 1,0)-cell or a (..., 1, l)-cell. In either case, Cj includes a cell of the form 
{(tt,x,^(M,x)) : u G D,gi{u) < x < hi{u)} for some semi-algebraic function 

Let c,d be real algebraic numbers such that gl < c < d < h^. Write c = cos{tp') and 
d = cos^tp) for 0 < -ip < Ip' < TT. By Lemmaj^the limits limt_>oo c) and limt_>oo C(e'’*, d) 
exist and are algebraic numbers in the interval [—1,1]. Let 6,8' G [ 0 , 7 r] be such that 
cos( 0 ) = limt^oo d) and cos( 6 i') = limt^oo C(e’’‘, c). 

By Corollary 3 we know that is either rational or transcendental. In particular 

we know that it is not equal to which is algebraic and irrational. Let us suppose that 
^ (the converse case is almost identical). Then there exists 6" with 9 < 6" < 9', 
such that 

9 <9"+^r^P<P'-iP) <9'. ( 14 ) 

bi 


Since 27 r, bi, 62 are linearly independent over Q it follows from Kronecker’s approximation 
theorem that {{bit,b2t) mod 2tt ■. t G K>o} is dense in [ 0 , 27 r)^ (see [T 31 Chapter 23 ]). Thus 
there is an increasing sequence ti < t2 < ■ ■ ■, with = tp mod 27 r for all n, such that 
b 2 tn mod 2tt converges to 9 ". Then, defining si < S2 < ... by s„ = -I- ^ , we have 
biSn = Ip' mod 27 r for all n and, by ( 141 , 


lim b2Sn = lim ^2^71 + ^{ip' — ’ 4 ’) = — ip) < 9 ' (mod 27 r) 

n—>-oo n—^00 Oi hi 


Let 77(f) = ^(e’'*,cos(&if)) — cos{b2t). Then for t G i?++ such that g{e'’') < cos(6if) < 

h(e'-‘), 

77(f) = 0 cos(62f) = ^(e*"*, cos(6if)) 

(e’’*,cos(6if),cos(62f)) G Cj 
tG Zj (by . 

Now lim„_>.oo ri{tn) = cos( 0 ) — cos( 9 ") > 0 and lim„_>.oo v(sn) < cos( 9 ') — cos( 9 ') = 0 . 
Moreover for n sufficiently large we have [f„, s„] C i?++. It follows that 77(f) has a zero in 
every interval [f„,s„] for n large enough. We conclude that Zj is unbounded. 
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4.2.3 Case III: < h^. 

This case is symmetric to Case II and we omit details. ◄ 
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I A I Hardness 

In this section we show that decidability of the Continuous Skolem Problem entails significant 
new effectiveness results in Diophantine approximation, thereby identifying a formidable 
mathematical obstacle to further progress in the unbounded case. 

Diophantine approximation is a branch of number theory concerned with approximating 
real numbers by rationals. A central role is played in this theory by the notion of continued 
fraction expansion, which allows to compute a sequence of rational approximations to a given 
real number that is optimal in a certain well-defined sense. For our purposes it suffices to 
note that the behaviour of the simple continued fraction expansion of a real number a is 
closely related to the (homogeneous Diophantine approximation) type of a, which is defined 
to be 



Let [ni,n 2 ,n 3 ,...] be the sequence of partial quotients in the simple continued fraction 
expansion of a. Then, writing K{a) := supj,>g n^, it is shown in |19l pp. 22-23] that L{a) = 0 
if and only if K{a) is infinite and otherwise 

K{a) < L{a)-^ < K{a)+2. 

It is well known that a real number algebraic number of degree two over the rationals 
has a simple continued fraction expansion that is ultimately periodic. In particular, such 
numbers have bounded partial quotients. But nothing is known about real algebraic numbers 
of degree three or more—no example is known with bounded partial quotients, nor with 
unbounded quotients. Guy m asks: 

Is there an algebraic number of degree greater than two whose simple continued fraction 
expansion has unbounded partial quotients? Does every such number have unbounded 
partial quotients? 

In other words, the question is whether there is a real algebraic number a of degree at least 
three such that L{a) is strictly positive, or whether L{a) = 0 for all such a. 

Recall that a real number x is computable if there is an algorithm which, given any 
rational e > 0 as input, returns a rational q such that 19 — 2:1 < e. The main result of 
this section is Theorem which shows that the existence of a decision procedure for the 
general Continuous Skolem Problem entails the computability of L{a) for all real algebraic 
numbers a. Now one possibility is that all such numbers L{a) are zero, and hence trivially 
computable. However the significance of Theorem]^ is that in order to prove the decidability 
of the Continuous Skolem Problem one would have to establish, one way or another, the 
computability of L{a) for every real algebraic number a. 

Fix positive real algebraic a, c and define the functions: 

fi{t) = e*(l — cos(t)) -I- t(l — cos(at)) — csin(at), 
f 2 {t) = e*(l — cos(t)) -I- t(l — cos(af)) -|- csin(at), 

f{t) = e*(l — cos(f))-I-t(l — cos(at)) — c| sin(af)| = min{/i(t),/ 2 (t)}. 

Then fi{t) and f 2 {t) are exponential polynomials. Moreover it is easy to check that the 
function f(t) has a zero in an interval of the form (T, oo) if and only if at least one of fi{t), 
f 2 {t) has a zero in (T, oo). 

We will first prove two lemmas which show a connection between the existence zeros of 
f{t) and the type L{a). We then will derive an algorithm to compute L{a) using an oracle 
for the Continuous Skolem Problem, thereby demonstrating our desired hardness result. 
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► Lemma 11. Fix real algebraic a,c and e € Q with a,c > 0 and e G (0,1). There exists 
an effective threshold T, dependent on a,c,e, such that if f{t) = 0 for some t>T, then 
L{a) < c/27r^(l — e). 

Proof. Suppose f{f) = 0 for some t > T. Define Si = t — 2Trm and 82 = at — 2'Kn, where 
m, n G N and 81,62 G [—tt, tt). Then we have 


n 

a - 

m 


|<^2 - aSi \ 

2 'Km 


We will show that for T chosen large enough, if /(t) = 0 for t > T then we can bound 1^21 
and |ai5i| separately from above and then apply the triangle inequality to bound 1^2 — o5i|, 
obtaining the desired upper bound on L(a). 

t — TT T — TT 

Define 0 < a < 1 by = (1 — £^). Since m > , for sufficiently large T we 

have 


27r 


27r 


t > 2 'K(m — 1) > 2 TTma . 


(15) 


Furthermore, since 0x^/2 < 1 — cos(a;) for |a:| sufficiently small, we may assume that T is 
large enough such that the following is valid for |a:| < tt: 


if 1 — cos(a;) < ctt/T then aa;^/2 < 1 — cos(a:). 


(16) 


We have the following chain of inequalities, where (*) follows from f{f) = 0 and e*(l — 
cos(t)) > 0: 

W c|sin(at)| c|sin((52)| c\S 2 \ 

1 — cos( 02 ) = 1 — cos(at) < -=- < -. 


It follows that 1 — cos((52) < cir/t and so by (161 we also have 

<^^2 /r \ 

— < 1 - cos((52) • 


Combining the upper and lower bounds on 1 — cos((52) and using (151, we have 
... 2 c 2 c c 

l<^2| < — < 


at 2'Kma'^ m7r(l — e^) 

We next seek an upper bound on |5i|. To this end, let T be large enough so that 


ce * < 


for t > T. 


V 2 aat / 

Then the following chain of inequalities holds: 

^ < 1 - cos(<5i) 

c| sin(52)| - t{l - 003 ( 82 )) 


(17) 


< ce 

< 


< 


f—y 

\ 2 aatJ 
! ce 
VdoTra^m 


{ valid for all |<5i| < tt } 

{ since f(f) = 0 } 
{ since | sin(52)|, | cos(<52)| < 1} 
{by (fig } 


{ by (151 } 
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It follows that 


|a(5i| < 


C£ 


7 rm(l — 

Finally, by the triangle inequality and the bounds on |a(5i| and |(52|, we have 

C+ C£ C 


n 

a - 

m 


|(52 - a(5i| ^ |(^2| + |a<^i| ^ 


27rm 27rm 27r^m^(l — 2TT'^m'^{l — e) ^ 

so the natural numbers n,m witness L(a) < c/27r^(l — e). 


► Lemma 12. Fix real algebraic a,c and £ G Q with a,c > 0 and e S (0,1). There exists 
an effective threshold M, dependent on a,c^£, such that if L{a) < c(l — e)/27r^ holds and is 
witnessed by natural numbers n,m with m > M, then f{t) = 0 for some t > 2 ttM. 

Proof. Select M large enough, so that c(l — £)/ttM < tt and 


if |a;| < c(l — £)/ttM, then (1 — £)|a:| < | sin(a;)|. (18) 

Suppose now that L{a) < c(l — e)/27r^, let this be witnessed by n, m G N with m> M and 
define t := 27rm. We will show that /(<) < 0. This suffices, because f{t) is continuous and 
moreover is positive for arbitrarily large times, so it must have a zero on [t,oo). 

Since L{a) < c(l — e)/27r^, we have jam — n| < c(l — e)/27r^m. Therefore, we can write 
at = 2TTam = 27rn + S for some 6 satisfying |(i| < c(l — £)/TTm < tt. We have 


fit) 

= { as cos(t) = 1 } 
t(l — cos((5)) — c| sin((5)| 

< { by (18) and 1 — cos(a;) < x'^/2 } 
TTmS^ — c(l — e)|(5| 

< { by |(5| < c(l — £)fTTm } 

0 . 


◄ 


The following corollary is immediate: 

► Lemma 13. Fix real algebraic a,c and £ G Q with a,c > 0 and e G (0,1). There exists 
an effective threshold T, dependent on a,c,£, such that if fit) 0 for all t>T, then either 
L(a) < c(l — e)/27r^ and this is witnessed by natural numbers n,m with m < T/2tt, or 
L(a) > c(l — e)/27r^. 

We now use the above lemmas to show the central result of this section: 


► Theorem 14. Fix a positive real algebraic number a. If the Continuous Skolem Problem 
is decidable then L{a) may be computed to within arbitrary precision. 


Proof. Suppose we know L{a) G [p, q\ for non-negative p,q GQ. Choose a real algebraic c 
with c > 0 and a rational e G (0,1) such that 


c(l — e) c 

2 tt '^ 27r^(l — e) 


< q- 
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Write A := c(l — e)/27r^ and B := c/27r^(l — e). Calculate the maximum of the thresholds T 
required by Lemmas 11 and 13 Check for all denominators m < T127: whether there exists 
a numerator n such that n, m witness L(a) < A. If so, then continue the approximation 
procedure recursively with confidence interval [p, A]. Otherwise, use the oracle for the 
Continuous Skolem Problem to determine whether at least one of fi{t), f 2 {t) has a zero on 
[T, oo). If this is the case, then f{t) also has a zero on [T, oo), so by Lemma [XT] L(a) < B 
and we continue the approximation recursively on the interval [p,B]. If not, then L(a) > A 
by Lemma 13 so we continue on the interval [^, 9 ]. Notice that in this procedure, one can 
choose c, s at each stage in such a way that the confidence interval shrinks by at least a fixed 
factor, whatever the outcome of the oracle invocations. It follows therefore that L{a) can be 
approximated to within arbitrary precision. ◄ 


We conclude by remarking that the exponential polynomials /i and /2 involved in the 
proof of Theorem |14| involve only two rationally linearly independent frequencies. Thus the 
theorem applies as soon as we have a decision procedure for exponential polynomials with 
two rationally linear independent frequencies. 


B I Proof of Lemma 

We divide Lemma into two separate results, which are proven below. 

► Proposition 15. There is a procedure that, given a semi-algebraic set D C and real 
algebraic numbers n,..., r„, returns an integer T such that {t > 0 : (t, ..., e'’"*) G D} 

either contains the interval (T, 00 ) or is disjoint from (T, 00 ). The procedure also decides 
which of these two eventualities is the case. 


Proof. Consider a non-zero polynomial P G IR[xo,..., x„] whose coefficients are real algebraic 
numbers. Then we can write P(t, ..., e’’"*) in the form 


for non-zero univariate polynomials Qi, ■ ■ ■ ,Qm with real algebraic coefficients and real 
algebraic numbers Pi > ... > It is clear that for t sufficiently large, P{t, ..., a’""*) 
has the same sign as the leading term Qi{t). The proposition easily follows. ◄ 

► Lemma 16. Let g : D ^ M. be a hounded semi-algebraic function with domain D C K". 
Let r = (ri,...,r„) be a tuple of real algebraic numbers and Ti an integer such that 
e'’* = ..., e’’"*) G D for all t > Ti. Then the limit g* = limt_,.oo 5 '(e’'*) exists, is an 

algebraic number, and there are effective constants T 2 ,e > 0 such that \g{e'^*) — g*\ < 
for all t > T 2 . 

Proof. Since g is semi-algebraic, there is a non-zero polynomial P with real algebraic 
coefficients such that P{x,g{x)) = 0 for all a; G I? (see 0 Proposition 2.86]). In particular, 
we have P(e’’*, 5 (e’'*)) = 0 for all t > Ti. Gathering terms, we can rewrite this equation in 
the form 

+ ... + = 0 


for non-zero univariate polynomials Qi, ■ ■ ■ ,Qm with real-algebraic coefficients and real 
algebraic numbers Pi > ... > pm- 

If m = 1 then for all t > Ti we have Qi(g(e^*)) = 0. Thus g{e‘^*) is equal to some of root 
of Q for all t > Ti. Then by Proposition 15 there exists T 2 such that g(e'^*) = g* for some 
fixed root g* of Qi and all t > T 2 . 
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If m > 1, since g is a bounded function, for all t > Ti we have 








(19) 


for some constant M. If Qi has degree d then (19) implies that the closest root of Qi to 
5 (e’’*) has distance at most Hence there exists a root g* of Qi and effective 

constants £,r 2 > 0 such that \g{e'^*) — g*\ < e“®* for all t > T 2 . ◄ 


I C I One Linearly Independent Frequency 


► Theorem 17. Let f{t) = J2j=i (-Pij(^) cos(wjt) + sin(wj)) he an exponential poly¬ 
nomial such that uji,... ,iOn o,re all integer multiples of single real algebraic number b. Then 
we can decide whether or not {t G IR>o : f{f) = 0} is bounded and, if bounded, we can 
compute an integer T such that {t G M>o : f{t) = 0} C [0,T]. 

Proof. Recall that for each integer n, both cos{nbf) and sm{nbt) can be written as polynomials 
in sin(&t) and cos{bt) with integer coefficients. It follows that we can write / in the form 

f{t) = Q{t, ..., e'’"*, cos( 6 t), sin( 6 t)), 


for some multivariate polynomial Q with real algebraic coefficients that is computable from 

/• 

Write K>o as the union of R+ := {t G K>o : sin(&t) > 0} and R- := {t G K>o : sin(&t) < 
0}. We show how to decide boundedness of {t G R+ : f{t) = 0}. The analogous case for 
R- follows mutatis mutandis. The idea is to partition {t G R+ : /(t) = 0} into components 
Zi,, Zm and to show for how to decide boundedness of each Zj. 

To this end define a semi-algebraic set 

E := {{u,x) G : 3?/ > 0 (cc^ = 1 A Q{u, x,y) = O)} . 

Then for all t G R+ we have f(t) = 0 if and only if (t, e“*, cos(bt)) G E. 

Consider a cell decomposition E = C\U .. . UCm, for cells C\,..., Cm C and define 


Zj = {t G R+ : (t, e'’*, cos(&t)) G Cj} 


for j = 1,..., m. Then {t G R+ : f{t) = 0} = U ... U Zm. 

We now focus on Zj for some fixed j G m} and show how to decide whether or 

not Zj is finite. Write Dj C for the projection of the corresponding cell Cj C 

onto the first n 3- 1 coordinates. Then Zj C {t G IR>q : (t, e’’*) G Dj}. 

First, suppose that {t G M>o : (Ce’’*) G Dj} is bounded. Since Dj is semi-algebraic, by 
Proposition |15| we can compute an upper bound T of this set. In this case t < T whenever 


t G Zj. 


On the other hand, suppose that {t G K : (t,e’’*) G Dj} is unbounded. Then, by 
Proposition 15 this set contains an unbounded interval (T, 00 ). We claim that in this case 
Zj must be unbounded. There are two cases according to the nature of the cell Cj. 


Case I: Q is a (..., 0)-cell. 

In this case there is a continuous semi-algebraic function g : Dj ^ M. such that Cj = 
{{u,g{u)) : u G Dj}. Then for t G i?+ H (T, 00 ), 

g(t, e'’*) = cos{bt) <;=^ (t,e’’*,cos( 6 t)) G Cj 

^ m = o. 
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M Figure 1 Intersection points of g{t, e’’*) and cos(fet) for t € -R+ (with shown in bold). 


In other words, / has a zero at each point t G i?+ n (T, oo) at which the graph of g(t, e’’*) 
intersects the graph of cos(bt). Since g is a continuous function with an unbounded domain 
that takes values in [—1,1], there are infinitely many such intersection points-see Figure]^ 

Case II: Cj is a (..., l)-cell. 

In this case Cj contains some (..., 0)-cell and so the argument in Case I shows that / has 
infinitely many zeros in Zj. ◄ 






